Spin Hall Effect in a Diffusive Rashba Two-dimensional Electron Gas 
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A nonequilibrium Green's functions approach to spin-Hall effect is developed in a diffusive two- 
dimensional electron system with Rashba spin-orbit interaction. In the presence of long-range 
impurities, the coupled quantum kinetic equations are solved analytically in the self-consistent 
Born approximation. It is shown that the intrinsic spin-Hall effect stems from the dc-field-induced 
perturbation of the density of states. In addition, there is an additional disorder-mediated process, 
which involves the transition of nonequilibrium electrons between two spin-orbit-coupled bands. It 
results in an additional collision-independent spin-Hall conductivity and leads to the vanishing of 
the total spin-Hall current even at nonzero temperature. 

PACS numbers: 72.10.-d, 72.25.Dc, 73.50.Bk 
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The discovery of spin-Hall effect, namely, the nonvan- 
ishing of spin current along the direction perpendicular 
to the dc electric field, has attracted much recent interest 
in the spin-dependent transport in two-dimensional (2D) 
semiconductors with the spin-orbit (SO) coupling. The 
early studies on this issue have been devoted to the spin- 
orbit interaction between electrons and impurities, 1 ' 2 
and therefore, the spin-Hall effect exhibits the extrin- 
sic character. More recently, the scattering-independent 
intrinsic spin-Hall effect, which entirely stems from 
the internal-field-induced spin-orbit coupling, has been 
predicted in the hole-doped semiconductors^ and two- 
dimensional semiconductors with Rashba 4 and Dressel- 
haus SO coupling. 5 

In clean two-dimensional electron systems (2DES) with 
Rashba SO interaction, Sinova et al. have shown that the 
spin-Hall conductivity, <t s h, has a universal value e/8ir 
at zero temperature^ Subsequently, a great deal of works 
has been focused on the effect of disorder on the spin-Hall 
current. When the two-dimensional systems are suffi- 
ciently dirty and the Anderson localization is dominant, 
the spin-Hall conductivity is found to be independent of 
the disorder £ In diffusion regime, the studies have been 
carried out when the presence of short-range disorders. 
It has been demonstrated in Ref. 

urn 

that the collision 

broadening can reduce the spin-Hall current and the bal- 
listic value of <7 s h can only be held for weak disorder. 
However, the further investigation revealed that the per- 
sistent spin-hall conductivity should be completely sup- 
pressed. This conclusion has been obtained by different 
methods, such as the Kubo formulajiSiii* 1 ^ the Keldysh 
formalism^ and spin-density method 1 ^ etc. It is com- 
monly accepted that such complete cancellation of spin- 
Hall current is not due to any symmetrjii^ and closely re- 
lates to the isotropy of the electron-impurity scattering^ 

However, for the realistic 2D systems, it is well known 
that the disorder collision being short-range is a crude 
approximation and the electron-impurity scattering is al- 
ways long-ranged. In Ref.0, by considering the effect 
of long-range disorders on the spin-Hall conductivity in 
2D electron systems with Rashba spin-orbit coupling, the 
nonvanishing of spin-Hall effect has been shown due to 



the anisotropy of the collision. This argument docs not 
take the full effect of collision anisotropy, however and 
hence a more careful investigation should be performed. 12 
In this letter, we construct a nonequilibrium Green's 
functions approach to the spin-Hall effect in the diffusive 
Rashba two-dimensional electrons. When the presence of 
the long-range impurity-electron scattering, the derived 
kinetic equations are solved analytically. We clarify that 
the dc-field-induced intrinsic spin-Hall current is com- 
pletely cancelled by the spin-Hall effect, stemming from 
the disorder-induced transition of the perturbative elec- 
trons, and in result, the total collision- independent spin- 
Hall effect also vanishes when the presence of long-range 
disorders. 

We consider a quasi-2D electron semiconductor in the 
x — y plane subjected to the Rashba SO interaction. The 
single-particle Hamiltonian of system can be written as 
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where a is the Rashba coupling constant, a are the Pauli 
matrices, m is the electron effective mass, and n is the 
unit vector perpendicular to the 2DES plane. By the 
local spinor unitary transformation [/(p) 11 



U(P) 

the Hamiltonian 
U+(p)h(p)U(p) 
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Q can be diagonalized, 
= diag(ei(p),e 2 (p)) with 



h'(p) 



p 2 /2m + (-iyap and fj, = 1, 2. 

The nonequilibrium Green's functions G r,< (p,ui) for 
the 2DES with Rashba coupling are defined in the gen- 
eral manner and become 2x2 matrices in spin spaced 
Obviously, in the helicity basis, where the single Hamil- 
tonian is diagonal, the noninteracting Green's functions 
are also diagonal 



Gq(p,w) = diag 



I 



1 



oj — E\ (p) + ie ' ijj — £2 (p) + ie J ' 



Gq(p, uj) — -2mF(w)ImG5(p,w) 



(3) 
(4) 
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Here np(w) is the Fermi function. 

In quasi-2D semiconductors, the electrons should ex- 
perience scattering by impurities. The previous stud- 
ies have only treated the short-range interaction between 
electrons and impurities and the relaxation is described 
by an isotropic parameter r. However, when the pres- 
ence of two spin-orbit-coupled bands, such a simplifica- 
tion loses many important relaxation times, namely, the 
longitudinal and transverse times. In this letter, we study 
more realistic long-range collision and distinguish these 
different scattering times carefully. 

In the spin basis, the electron-impurity scattering 
process is described through the long-range potential 



V(p — k). By transforming to the helicity basis it be- 
comes f (p, k) = U+(p)V{p - k)L>(k). To deal with the 
electron-impurity scattering, we restrict ourselves to the 
self-consistent Born approximation, which corresponds to 
the vertex corrections in the scheme of Kubo formalism. 
Hence, the self-energies can take the form 



(P,W) 



k 



f(p,k)G r '<(k,c l ;)f+(p,k) ) (5) 



with the impurity density rii 
Eq. JSJ gives 



Putting the matrix U into 
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E r '<(p, T, t) = l -n, \V(P - k)| 2 { ai G r '< + a 2 a x G r <<a x + ia 3 [a x , G r ><}} 
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Here at(i — 1,2,3) are the factors associated with the 
directions of the momenta, ai = 1 + cos(0 p — 0k) , a i — 
1 - cos(0 p - k ), a 3 = sin(0 p - k ). 

In order to investigate the transport properties of 2D 
system driven by a dc electric field E along x axis, the 
dominant task is to carry out the nonequilibrium less 
Green's functions G < . The measurable quantities such 
as current and spin-current, are closely related to them. 
It is well known that the electric current operator in the 
spin basis reads J = ^ p with the one-particle 

matrix current operator j being j = e \j- + a(n x er)] . 
The current operator of the z-direction spin can be de- 
fined through the electric current operator: P = (J<Ji + 
<7iJ)/4e. 17 We find, in the helicity basis, the operator of 
z-direction-spin current along y axis is nondiagonal. By 
taking the statistical average, the Hall spin-current can 
be computed via 



E Py 
2m 



[Pi2 (p) +P2l(p)] 



(7) 



with the distribution functions Pap{p) — 
-*/^G< /3 (p,w)(a,/3 = 1,2). It can be clearly 
seen that the presence of spin-Hall effect is due to 
the nonvanishing of the p±2 and p%\, viz. interband 
polarizations. 

We follow the standard procedures described in Ref.lTsI 
to derive the kinetic equations. We first construct the 
equations for the less Green's functions in the spin basis, 
and then modify it by the local unitary transformation. 
Finally, for steady and homogeneous systems, we arrive 
at the kinetic equations in the helicity basis 
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with B K = [G < ,cr 2 ]_. In this formula, the arguments of 
the Green's functions (p, to) are dropped for shortness. 
In deriving the equations, the scalar potential gauge is 
used and only the lowest order of gradient expansion is 
taken into account. Noted that the second term in the 
bracket of the left hand of Eq. (JSJ arises from the fact 
that the transformation is local. As we will show below, 
it will lead to the intrinsic spin-Hall conductivity. At the 
same time, when the electron-impurity collision is cho- 
sen to be short-ranged, these kinetic equations reduce to 
the formula obtained in Ref.^1 by means of the Keldysh 
formalism. 

We are only concerned with the linear response of the 
system to dc fields. Hence, the kinetic equations can be 
linearized with respect to the external electric fields. We 
rewrite the less Green's functions as G < = Gq + Gf, 
with Gq and Gf being in the zero- and first-order of E, 
respectively. As a result, the kinetic equations for the 
less Green's functions Gf have the forms 



a P C! + ieE ■ V P G< - ^eE • V p p A) = 



(S r G< - G<S a - G r S< + £<G a )i, 



where the matrices C\ and D are 
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(9) 



(10) 



(11) 



The subscript 1 in the right-hand side of Eq. © stands 
for keeping this term in the first-order of E. 

From Eq. 10 we can see that two driving terms enter 
the kinetic equations: ieE • VpG^ and — eE ■ V p 0p-Do/2. 
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Hence we can assume that the solutions of these equa- 
tions are formed from two terms and Gj 2 . In fact, 
when the disorder becomes short-ranged, these two parts 
of solutions connect with the terms in Kubo formalism. 
The first part of solutions stems from the driving term 
completely and corresponds to the bubble diagram in the 
scheme of Kubo formulaii2iiiii2 At the same time, the an- 
other part of solutions, relating to the electron-impurity 
scattering, should be understood as the vertex correction. 
In this way, the distribution functions in the first order of 
dc field, pi(t), can be written as pi(t) — Pi (t) + p{\t) 
and, correspondingly, the spin-Hall effect comes from two 
different interband-polarization processes. 

The less Green's functions G^ , associated with the 
matrix Do, relate to the function n-p(uj). It is a non- 
diagonal matrix with same elements. In fact, Gj 
originate from the nondiagonal elements of linear re- 
sponse retarded Green's functions, (Gi )r2 = (61)21 = 
— 2mFlm(G£)i2. The latter one, (61)12, obeys the Dyson 
equation 

2op(GI)ia + ^eSsin0 p [(GS) 11 - (G r Q ) 22 } = 
zp 

(GI) ia [(E5) u - (E5) M ] - (EDi2[(G5) u - (G£) 22 ],(12) 
and can be obtained analytically as 

eE 

{G\)i2 = sin^ p [(GS)u - (G5) M ]. (13) 

It is well known that the linear response can not dis- 
turb the retarded Green's functions in semiconductors 
without spin-orbit interaction, — However, when the spin- 
orbit coupling is introduced, even weak dc field can 
produce the transition between the spin-orbit-coupled 
bands, resulting in interband polarization. The factor 
np(oj) in the less Green's functions Gi indicates that 
all electrons join in this process. 

Substituting the Eq. I|13|) into and neglecting the 
broadening of the noninteracting retarded Green's func- 



tions, we find the contribution from the first polarization 
process to the spin-Hall conductivity 

a sH = Ta / dpn>F(£i(p) - n) - n F (e 2 (p) - p), 

where p is the chemical potential. At zero tempera- 
ture, a^2 approaches the ballistic value e/8n. Note that 
this formula agrees with that obtained in the previous 
studies! 7 ! 1 ^ 12 

The remaining part of Eq. @ describes the transport 
process and is associated with the function dnF(ui)/duj. 
To solve it, we employ the two-band generalized 
Kadanoff-Baym ansatz (GKBA). This ansatz is widely 
used in the quantum kinetics of semiconductors driven by 
ac field o 20 ' 21 and accurate enough to yield the quantita- 
tive agreements with experiments^ To study the linear 
response of 2DEG with SO coupling, we only need to con- 
sider the GKBA in the first order of dc field strength E. 
To further simplify the treatment, we ignore the broad- 
ening of Gq a induced by electron-impurity scattering. It 
is well known, this approximation combined with the 
gradient approximation are referred as the Boltzmann 
limit and valid when the spatial and temporal variations 
change slowly^ They are effective for modelling the qua- 
siclassical optical and transport processes in semiconduc- 
tors. 

By taking into account the time-independence of p(t) 
in the steady-state transport, the second part of solution 

of Eq. (O can be obtained analytically. Retaining the 

(2) 

results to O(rij), we find that the functions (p\ )i 2 can 
be expressed as 

(/0i 2) )i 2 (p) = C(p)eE$m(f)p + i£(p)eEsm(l) p , (15) 

where £(p) and £(p) are real functions. Only the real part 

of (p\ )i 2 needs to be calculated, since the contribution 
of the imaginary part to the spin Hall current vanishes 
due to the symmetry relation (p\ )i 2 = (p x ) 2 i. 
The elementary calculation yields 



CO) = |Vjsnp(i7) \ E = S2 {p)(a + p/m) - W E n F (E)\ B=e ^ p) (-a + p/m) 



Aap 

H — — [02 (p) + ©1 (p+ 2ma)] (a+p/m) — [0i(p) + 2 (p- 2ma)] (-a+p/m) } . (16) 

T221 T212 



I 

The functions 0i(p) and © 2 (p), coupled by the equa- connect to the diagonal distribution functions 
tions, 



e p (p+(-l)"2ma) 

-V En F (E)\ E=e ^ p ) ~ 



OiW(p) = [(-l)"a+p/m]e^(p)e£7coB^p, (18) 



(17) with /i = 1, 2 and p, = 3 — fx. In these equations, several 
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different relaxation times, defined by 

— = 2nm \V(P - k)| 2 A VI ,(^ k _ p ,p, k), (19) 
T%lLV k 

emerge due to the long-range potential. Here the 
are A lfJ , v (<p,p,k) = \ sin 2 <M(e w - ^k), ^2^((f>,P,k) = 
|(1 - cos^)(5(e w - £„ fc ) and A 3fiv ((f),p,k) = |cos^(l - 
cos^)<y(£pp - e„fc). 

The contribution of second intcrband polarization pro- 
cess to the spin-Hall conductivity can be evaluated by 

'^/"^CGOdp. (20) 

It should be noted that the spin-Hall conductivity com- 
ing from the terms in the second line of Eq. (|16(l is zero 
at zero temperature. This is straightforwardly educed 
from the fact that zero-temperature difference between 
the Fermi momenta for the two spin-orbit-coupled bands 
is 2ma and the factor \7 E n F (E) reduces to the delta 
function. Omitting these terms at nonzero tempera- 
ture is also supported by numerical estimation. At tem- 
perature lower than 4.2 K, considering long-range col- 
lision between the remote impurities and electrons in 
GaAs/AlGaAs heterojunctionsr 2 " we find the magnitude 
of contributions of these terms to the spin-Hall conduc- 
tivity more than five orders smaller than that from the 

(2) 

other terms. Hence, the spin-Hall conductivity a sI ^ can 
be simplified as 

= Wnma J ?dp ^ eUf ^ U=^(p)(» +P/ m ) 
-V E n F (E)\ E=ei{p) (-a+p/m)}. (21) 

Finally, we obtain — — c^y. 

From the Eqs. (|16fl and (|17|) we can see that the spin- 

(2) 

Hall conductivity a sI ^ appears from the transport pro- 



cess. In the dilute impurity limit, the diagonal distribu- 
tion functions of nonequilibrium electrons driven by dc 
fields are of order of n^ 1 . These perturbed electrons in 
one band can transit to another band, resulting in in- 
tcrband polarization. The latter process is proportional 
to rii and, therefore, the nondiagonal distribution func- 
tions in the lowest order of rij appear. Hence, it is ob- 
vious that, here, the disorder plays an intermediate role. 
Such disorder-mediated process, in contrast to the one 
produced by all electrons, contributes from the electrons 
near the Fermi surfaces. 

Note that in our treatment the influence of broadening 
of Green's functions on the spin-Hall effect remains un- 
touched. This implies that our results are useful for the 
relative clean samples. At the same time, the derived ki- 
netic equations in this letter can serve as the start point 
for further studying the effect of collision broadening, 
which is expected to affect on the linear response of the 
retarded and less Green's functions in different fashions. 
Hence, the spin-Hall effect may be observed in dirty sam- 
ples. 

In conclusion, we have presented a nonequilibrium 
Green's functions approach to the spin-Hall effect in 
Rashba two-dimensional electron systems. By taking 
into account the long-range scattering between impuri- 
ties and electrons, the intrinsic spin-Hall conductivity 
has been recovered. It has also been illustrated that the 
interband transition of nonequilibrium electrons, partic- 
ipating in longitudinal transport, results in a disorder- 
mediated spin-Hall effect. As the result, the total spin- 
Hall effect disappears. 
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